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Abstract. Accretion rates from disks around pre-main sequence stars are of importance for our understanding of planetary 
formation and disk evolution. We provide in this paper estimates of the mass accretion rates in the disks around a large sample 
of Herbig Ae stars. We obtained medium resolution 2 fim spectra and used the results to compute values of M acc from the 
measured luminosity of the Br r emission line, using a well established correlation between L(Br y ) and the accretion luminosity 
Lace- We find that 80% of the stars, all of which have evidence of an associated circumstellar disk, are accreting matter, with 
rates 3 x 10~ 9 < M acc < 10~ 6 M Q /yr; for 7 objects, 6 of which are located on the ZAMS in the HR diagram, we do not detect any 
line emission. Few HAe stars (25%) have M acc > 10~ 7 M Q /yr. In most HAe stars the accretion rate is sufficiently low that the 
gas in the inner disk, inside the dust evaporation radius, is optically thin and does not prevent the formation of a puffed-up rim, 
where dust is directly exposed to the stellar radiation. When compared to the M acc values found for lower-mass stars in the star 
forming regions Taurus and Ophiuchus, HAe stars have on average higher accretion rates than solar-mass stars; however, there 
is a lack of very strong accretors among them, probably due to the fact that they are on average older. 



Oh' 

6 

CO 



X 



1. Introduction 

Pre-main sequence stars accrete matter from their circumstellar 
disks for a large fraction of their life, until the disk is accreted 
entirely onto the star or is dissipated by a number of potentially 
competing processes, such as dynamical perturbations by close 
companions, photoevaporation by the host or nearby stars, or 
planet formation. Although the amount of matter accreted is 
small compared to the mass of the central object, accretion 
of matter through the disk has important consequences on the 
physics of the disk itself and on the properties of the planetary 
systems that the disk may form. 

The mass accretion rate through the disk is reasonably well 
known for pre-mains sequence stars with masses from about 
one solar mass (T Tauri stars or TTS), down to brown dwarfs 
(BDs). Results can be found, e.g., in Gullbring et al. 1998 
Muzerolle et al. 119951120051120351 . W hite an d Ghez COOfl . 
Natta et al. 120041 12006L Mohanty et al. ( 12005b . 

It has been much more difficult to determine M acc in stars 
more massive than about 1 M . Recently Calvet et al. @004) 
obtained the accretion rate in nine young pre-main sequence 
stars with masses in the interval 1.5-3.7 M (intermediate mass 
TTS or IMTTS) using HST and optical spectra. They found an 
average accretion rate of ~ 3 x 10~ 8 M /yr, with a spread of an 
order of magnitude at least. Muzerolle et al. (|2004) derived a 



similar value of M acc for the ~ 2 M star UX Ori by fitting its 
Hq, profile with magnetospheric accretion models. 

We report in this paper measurements of the accretion rate 
in a large sample of Herbig AeBe objects (intermediate mass 
pre-main sequence stars), mostly of spectral type A (HAe). 
We derive the accretion rate from the luminosity of the Br r 
line seen in emission in pre-main sequence accreting objects. 
Calvet et al. (2004) have shown that the correlation between 
the hydrogen IR line luminosity and the accretion luminosity 
found in TTS and BDs (Muzerolle et al. [T9951 Natta et al.|2004 1 
extends to their sample of IMTTS (see also van den Ancker 
2004). Establishing the luminosity of the hydrogen near-IR re- 
combination lines as proxy of the accretion luminosity has been 
extremely important, since it has made possible to derive accre- 
tion rates with reasonably good accuracy in objects where mea- 
suring the UV veiling is challenging. This is the case, for ex- 
ample, of very young star-forming regions, such as Ophiuchus, 
where the high extinction is the limiting factor (e.g., Natta et 
al. 2006 ), but also of the Herbig AeBe stars, where it is difficult 
to separate the UV excess emission from the hot photospheric 
continuum. 

The main purpose of this work is to estimate M acc in ob- 
jects with disks for which we are currently obtaining near-IR 
interferometric data with the ESO VLTI (AMBER GTO 1 ; see 
also Malbet et al. 120051 Tatulli et al. 120061 . Knowledge of M acc 



Send offprint requests to: A. Natta, natta@arcetri.astro.it 



1 http://www.eso.org/observing/proposals/gto/amber/ 



2 



R. Garcia Lopez, A. Natta, L. Testi and E. Habart: Accretion Rates in Herbig Ae stars 



is crucial for understanding the disk geometry, in particular the 
nature of the disk inner rim (e.g., Natta et al. 2001 Dullemond 
et al. 120011 Isella & Natta 2005), which can only form if the 
accretion rate is lower than a critical value of the order of 10~ 7 
M /yr (Muzerolle et al. 2004 1. For low M acc , the inner gaseous 
disk is optically thin and does not prevent the direct stellar ir- 
radiation of the inner dust edge of the disk, thus giving raise 
to the pufFed-up rim. Secondly, we are interested in comparing 
the results to the accretion properties of the lower-mass ob- 
jects studied in Taurus and Ophiuchus. In both regions it has 
been found that the mass accretion rate M acc is a strong func- 
tion of the mass of the central object, increasing roughly as 
M* 1 - 8-2 - (e.g., Natta et al. 2006 and references therein; but see 
also Clarke & Pringle 2006 1. The causes of this behaviour are 
not understood, and are matter of debate (Muzerolle et al. 2005 
Natta et al. 120061 Alexander & Armitage 120061 Dullemond et 
al. 20061. The HAe stars are scattered in different regions of 
the sky and are, on average, older than objects in Ophiuchus 
and Taurus; they may provide interesting information on the 
dependence of the disk accretion on age and environment. 

The outline of this paper is as follows. The observations 
and data reduction are described in §2; method and results are 
presented in §3 and discussed in §4. Conclusions follow in §5. 

2. Observations and data reduction 

Medium resolution (R~9000) spectra of all stars in our sam- 
ple were obtained with the ISAAC spectrograph at the Antu 
8.2m VLT unit telescope in service mode in spring 2004. The 
spectral region covered by our spectra included the Bry and 
the H2 (1-0)S(1) lines, which was not detected in any of the 
objects of the sample. All spectra were obtained in moderately 
poor conditions (seeing ~ 1 .5 arcsec and non photometric trans- 
parency), as these were adequate for the scope of the project. 
As a consequence, we did not attempt to obtain an accurate flux 
calibration of the spectra. 

The raw data were wavelength calibrated using daytime 
arcs and rectified using data from the ISAAC calibration plan. 
After aligning and coadding the dithered integrations on each 
target, the spectra were extracted and telluric absorption and 
instrumental throughput were corrected using observations of 
telluric standards obtained by the observatory staff during the 
same nights as the science observations. The details of the 
ISAAC calibration plan, the reduction process and the instru- 
ment characteristics can be found on the ESO web pages 2 . 

The portion of the spectra around Br y are shown and briefly 
discussed in the appendix. In each spectrum we measured the 
equivalent width of the Bry line; the measurements are given 
in Tables 1 and 2, together with the lcr uncertainties. 

Tables 1 and 2 report also the star name and properties. 
Spectral type and distance are compiled from the literature 
(e.g., van den Ancker et al.[T997|[T998l Testi e t al.[T998l Meeus 
et al. l200Tl Rodgers I200T1 Hernandez et al. 120051 and refer- 
ences therein). Extinction and luminosity are computed from 
the V, (B-V) magnitudes as described in Testi et al. J1998i : 
stellar masses are derived by comparing the location of each 

2 http://www.eso.org/instruments/isaac 



star on the HR diagram to the evolutionary tracks of Palla and 
Stahler (JT993}. 

3. Method and results 

We derive the accretion luminosity L acc of each star from the 
luminosity of Br y , following the relation derived by Calvet et 
al. l|2004) for a sample of objects ranging in mass from about 4 
M Q to BDs: 

log L acc /Lo = 0.9 x (log(L(Br r )/Le) + 4) - 0.7 (1) 

where L(Br y ) is the luminosity of the Br y line emitted by the 
circumstellar gas. Note that Eq.(l) is an empirical relation, 
where L acc is measured from the observed veiling, i.e., inde- 
pendently of any assumption on the origin of the Br y emission. 
It has been derived from objects with L(Br y ) between ~ 3x 10~ 6 
and ~ 3 x 10~ 3 L ; the scatter of L acc for any given L(Br r ) is 
typically a factor +3 (see also van den Ancker 20041. 

In order to compute L(Br y ), we need to extract the equiv- 
alent width of the circumstellar component of the line from 
the observed one. In HAe stars, the observed flux is the sum 
of the stellar photospheric emission, which for A stars shows 
broad and rather strong absorption in the hydrogen lines, of 
the emission of the circumstellar gas and of the emission of 
the disk. Following Rodgers (2001), we assume that the star 
dominates the observed V emission. The equivalent width of 
the photospheric component, corrected for the K-band veiling, 
is computed for each object with spectral type earlier than A8 
using the template spectra of Rodgers (2001 1. For objects of 
later spectral type, for which we did not have spectra of tem- 
plate stars of suitable resolution, we compute W(Br y ) f ,y e from 
the expression: 

W(Br y ) circ = W(Br y ) obs - W(Br 7 ) phot 10-° AAK (2) 

where W(Br y ) OJ i,.s is the observed equivalent width, W(Br y ) p / iof 
the equivalent width of the photospheric Br y line of a template 
with the same spectral type (see Table |5J values from Rodgers 
2001 and Wallace & Hinkle 1997 1 and AK the disk continuum 
emission, computed subtracting from the observed K magni- 
tude the photospheric one. All the relevant parameters used 
in the derivation are given in Tabled K magnitudes are from 
2MASS 3 ; note that none of the quantities used in eq.(2) has 
been measured simultaneously. 

Of the 32 targets, 4 are very luminous objects of spectral 
type early B. Three of them have Br y in emission, one in ab- 
sorption. Table [3] refers in Column 3 a lower limit to the line 
luminosity, derived from the observed Br y equivalent width, 
i.e., without correcting for the photospheric absorption. For the 
three stars where it is seen in emission L(Br y )>0.1 - 0.2 L . 
The validity of Eq.(l) for these very luminous objects is very 

3 This publication makes use of data products from the Two 
Micron All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 
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Table 2. Herbig Be properties - Column 1 : star name; Column 
2: spectral type; Column 3: observed Br r equivalent width; 
Column 4: Br y luminosity. 



(1) 


(2) 


(3) 


(4) 


name 


ST 


Br y (obs) 


LogL(Br y ) 






(A) 


(L G ) 


V921 Sco 


BO 


-17.3±0.25 


>-0.7 


MWC 300 


B0-B1 


-10.9±0.25 


>-0.8 


MWC 166 


BO IV 


2.5±0.2 




MWC 297 


B1.5 


-16.0+0.2 


>-0.6 



doubtful, and we will not consider them further in the following 
of this paper. 

After subtraction of the photospheric component, 24 of the 
28 HAe stars have circumstellar Br y emission; for these, we 
compute L(Br y ) from W(Br y )„> c , the continuum flux derived 
from the extinction-corrected K magnitude and the distance of 
the star. L acc is then obtained from Eq.(l). For the remaining 4, 
where no circumstellar emission is detected, we estimate that 
W(Br r )„> c < 2 A, and compute upper limits to L acc accordingly. 
For HD 98922 and HD 144432, L acc refers to the quoted dis- 
tance lower limit, and is therefore also a lower limit to the true 
value. The results are given in Tableland shown in Fig.^ as a 
function of the stellar luminosity. Note that Eq.(l) has been de- 
rived for objects with L(Br y )<3 x 10 3 L ; 2 objects (HD 98922 
and VV Ser) have larger L(Br y ), and the use of Eq.(l) should 
be viewed with some caution. 

The mass accretion rate follows from the relation M acc - 
LaccR-kliPM*), and is shown in Tabled Column (13). 

The uncertainties on L acc and M acc are dominated in most 
cases by the uncertainties deriving from the subtraction of the 
veiled photospheric absorption component from the observed 
spectrum. The error bars shown in Fig^correspond to an un- 
certainty on the resulting W(Br y ) e „ e of ~ +2 A, which we con- 
sider realistic in the majority of cases. 

For any given object, variability of lines and continuum can 
also be large. There are no systematic studies of time variabil- 
ity of the Br y equivalent width in our sample stars, but one 
can compare our measurements with the Rodgers (2001 1 re- 
sults from lower resolution spectra (R~ 500 - 2000) for the 7 
stars we have in common. The differences are within the errors 
for 3 of them, but larger than 2 A in 2, one being R CrA. This 
confirms the well known fact that the accretion rate in pre-main 
sequence stars varies and that one can only estimate the "typ- 
ical" value of any individual object if the results of long term 
time monitoring are available. However, if one is interested in 
the properties of a large sample of stars, as in this work, one 
can safely assume that their time variability will increase the 
spread of the determinations, but not bias the average values. 

4. Discussion 

Fig-fflplots L flcc as a function of L* for the HAe observed in 
this paper. We have added six HAe stars (listed at the bottom 
of Table 1) taken from the Rodgers ( 2001 1 thesis, for a total of 
36 objects. Typically, the HAe stars have L acc roughly between 
10% and 1% L*. There are very few objects with L flee >0.1 L+, 
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Lg L + /L 

Fig. 1. Accretion luminosity as a function of L* for our sample 
(filled squares) and the Rodgers ( 2001 1 objects (open squares); 
arrows show upper limits. For comparison, we show also the 
IMTTS L acc measured by Calvet et al. 120041 from veiling 
(crosses). 

Table 3. Br y equivalent width of template stars - Column 1: 
spectral type; Column 2: adopted Br r equivalent widht. 



(1) 


(2) 


ST 


W(Br y ) 




(A) 


B8 


11.4 


AO 


9.7 


A2 


13.4 


A8 


12.0 


F3 


6.6 


F6 


5.5 



contrary to what happens for lower mass objects in Taurus and 
Ophiuchus (see Natta et al. 120061 Calvet et al. l2004t . 

Fig. |2 top Panel, shows the M acc distribution of the HAe 
stars; the values of M acc cover the interval ~ 10~ 6 -<3x 10~ 8 
M /yr, but few stars (25%) have M acc > 10~ 7 M /yr; 36% have 
M acc < 10~ 8 M /yr. The median value is M flcc = 3.5 x 10~ 8 
M /yr. These results agree with those obtained by Muzerolle 
et al. (2004), who found that intermediate-mass objects have 
typical H a profiles and/or Balmer discontinuity roughly consis- 
tent with the predictions of magnetospheric accretion models 
with accretion rates <10~ 7 M /yr. They studied in detail only 
one object, UX Ori, for which they found a value M acc ~ 10~ 8 
M /yr, comparable with the value (~ 6x 10~ 8 M /yr) we derive 
from the low resolution spectrum of Rodgers ( 2001 1. 
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Table 1. Herbig Ae properties - Column 1 : star name; Column 2: spectral type; Column 3: V magnitude; Column 4: K magnitude; 
Column 5: V extinction; Column 6: distance; Column 6: effective temperature; Column 8: stellar luminosity; Column 9: stellar 
mass; Column 10: observed Br y equivalent width; Column 11: equivalent width of the circumstellar Br y ; Column 12: accretion 
luminosity; Column 13: mass accretion rate. 
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Note: The last 9 lines are objects from Rodgers (2001 1. 



Our results support the conclusion of Muzerolle et 
al. ([2004 ) that in most HAe stars the accretion rate is suffi- 
ciently low that the gas in the inner disk, inside the dust evapo- 
ration radius, is optically thin, allowing the direct irradiation of 
the inner dust edge (the rim) by the star. One of the conditions 
for the formation of a puffed-up rim (Natta et al. I2001i seems 
therefore generally verified. 

Fig.[3]plots M acc as a function of M+ for the HAe stars, to- 
gether with objects in Ophiuchus (Natta et al. 2006 l, Taurus 
(Muzerolle et al. 120051 and references therein; Mohanty et 
al.|2)05 1 and IMTTS in Orion and Taurus (Calvet et al. l20l)4l 
Fig- EJcompares the M acc distribution of the HAe stars to that 
of objects in Ophiuchus in two mass ranges, as labelled. One 



can see that the median M acc for the HAe stars (~ 3 x 10~ 8 
M /yr) is about a factor ten higher than for Ophiuchus stars 
with 0.3 < M* < 1 M , and a factor 100 larger than in the very 
low mass objects with M*< 0.3 M . This is roughly consistent 
with the predictions of the M flec oc M* 1 8 relation derived for 
sub-solar mass objects. A close inspection of Fig. O however, 
shows that the HAe stars lack very strong accretors, as already 
noted in discussing Fig.^ If the upper envelope of the distribu- 
tions observed in Taurus and Ophiuchus extended to HAe stars, 
we would expect a much larger fraction of objects accreting at 
rates higher than 10~ 7 M Q /yr. 

It is not easy, at this stage, to understand why this happens. 
The HAe stars form a sample which is not complete nor ho- 
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Fig. 2. Distribution of the mass accretion rate in HAe stars (top 
Panel). The shadowed boxes show actual measurements, the 
empty boxes (thick histogram) the distribution of the upper lim- 
its. The mid and bottom panels show the distribution of M acc 
for Class II objects in Ophiuchus with mass 0.3 < M*/M Q < 1 
and MJM Q < 0.3, respectively (data from Natta et al. 2006 1. 

mogeneous. The objects are located in a variety of star forming 
regions, and, in some cases, are isolated in the sky. We think 
that the most likely interpretation of their behaviour is "ag- 
ing", since the accretion rate in viscous disks with constant a 
parameter is expected to deacrease with time roughly as r 15 
(Hartmann et al. 1998 1. HAe stars are, on average, older than 
Ophiuchus objects, which cluster at ages ~ 0.5 - 1 Myr, re- 
spectively (Palla & Stahler 2000 1. Their nominal ages, obtained 
by comparing their location on the HR diagram with Palla & 
Stahler (1999 1 evolutionary tracks, range from ~ 1 to > 10 
Myr; several stars are close to the ZAMS, and it is in fact inter- 
esting that 6/7 stars with non-detected Br r emission are ZAMS 
objects, the exception being HD 35187, with an estimated age 
~ 8 Myr. However, one should not overinterpret these findings, 
as ages of intermediate-mass stars derived from their location 
on the HR diagram are very uncertain. 

It is possible that the steep dependence of M acc on M* does 
not extend to objects more massive than, say, ~ 1 M , reflecting 
a change in the physical processes that control disk evolution 
and/or formation, and which have been discussed in a number 
of recent papers (e.g., Muzerolle et al. 2003 Natta et al. 2006 
Padoan et al. 120041 Alexander & Armitage 120061 Dullemond 
et al. 2006 1. Unfortunately, given the lack of completeness and 
the selection biases that affect it, it is difficult to establish the 
relevance of the HAe sample for this ungoing discussion. It is 
worth to remind, in this context, the well known bias due to the 
spectral type selection, which introduces an age-mass correla- 



Fig. 3. Mass accretion rate as a function of M* for HAe stars 
(filled squares and upper limits), Ophiuchus Class II objects 
(open dots and upper limits; Natta et al. 2006), Taurus CTTS 
and BDs (crosses; Muzerolle et al. 120051 Mohanty et al. 120051 
Calvet et al. f2004| and references therein), IMTTS in Orion 
(plus signs, Calvet et al.|T004). 

The dashed line shows the best-fit to the Ophiuchus data (M acc oc 
M* LS ); the dotted lines are ±1 dex. 

tion in the sample (e.g., van Boekel et al. 2005 1 and automati- 
cally excludes young, relatively low mass (< 2.5 M ) objects. 

5. Conclusions 

We have reported in this paper measurements of the disk mass 
accretion rates in a sample of 36 Herbig Ae stars. The values 
of Macc are derived from the measured luminosity of the Br y 
emission line, using the correlation between L(Br y ) and the ac- 
cretion luminosity L acc , established by Muzerolle et al. ( 1998 1 
and Calvet et al. ( 2004 1. Note that it is an empirical correlation, 
which makes no assumptions on the origin of Br y . 

We find that 80% of the stars, all of which have evidence 
of an associated circumstellar disk, are accreting matter, with 
rates 3 x 10~ 9 <M acc < 10~ 6 M /yr; for 7 objects, 6 of which are 
located on the ZAMS in the HR diagram, we do not detect any 
line emission. Few HAe stars (25%) have M acc > 10" 7 M Q /yr. 
These results support the conclusion of Muzerolle et al. (2004i 
that in most HAe stars the accretion rate is sufficiently low that 
the gas in the inner disk, inside the dust evaporation radius, 
does not prevent the formation of the puffed-up rim, proposed 
by Natta et al. (200D and Dullemond et al. (200T). 

When compared to the M acc values found for lower-mass 
stars in the star forming regions Taurus and Ophiuchus, HAe 
stars have on average higher accretion rates; however, there is 
a lack of very strong accretors among them, as shown by the 
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very few objects with L flcc > 10% L*. It seems likely that this is 
because on average HAe stars are older than TTS in Ophiuchus 
or Taurus. However, ages of intermediate-mass stars derived 
from the HR diagram location are very uncertain; this, and the 
various selection biases which plague the HAe sample, make 
the comparison very difficult. 

When we consider individual stars, the largest uncertain- 
ties on Mace come from line variability, as in TTS, but also 
from the presence of a photospheric absorption component of 
Br r in early-type stars, which needs to be subtracted from the 
observed one, and from the K-band excess emission. A better 
characterization of the photospheric profiles of template stars 
would be very valuable. 
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Appendix A: Br, Spectra 

The observed spectra, normalized to the continuum level and 
shifted for an easier display, are shown in Fig. IA.1I They are 
ordered roughly according to the observed W(Br r ) ois , from 
strong emission to strong absorption. In HAe, the strongest 
peak emission (HD 97048) is about 44% of the continuum; the 
deepest absorption (HR 4796) is 15% of the continuum. 

The resolution of the spectra (~ 30 km s _1 ) allows us 
to resolve the lines in all objects. The observed line shape 
results from the superposition of the circumstellar emission, 
the photospheric absorption and the veiling K-band contin- 
uum, very likely from the disk itself. One can see the signa- 
ture of the broad photospheric absorption component in sev- 
eral objects, where the K-band veiling is moderate or absent 
(e.g., HD 149914, HD 97300, HR 4796). In some objects both 
the circumstellar emission and the photospheric absorption are 
clearly visible (e.g., HD 141569 and 51 Oph); in both objects 
the emission has a double-peaked profile with peak-to-peak 
separation ~ 300 km s . 

In general, the emission lines are very broad, with half 
width at 10% of peak intensity (10%HW) of 200-300 km s" 1 . 
VV Ser has very broad wings, which extends to + 450 km 
s _1 from the line center. As a reference, the escape velocity 
of a typical AO star (M*=2.3 M , R=4 R Q ) is 330 km s" 1 ; the 
10%HW of Br r is often close to this value. 

In comparison, the three early B objects with emission have 
very strong lines. They are relatively narrow (10%HW <120 
km s ) when compared to the estimated escape velocity, which 
are of ~ 200 km s _1 . 

Clearly, the study of the profiles of the near-IR hydro- 
gen emission lines in these intermediate-mass objects promises 



very interesting results, and should be pursued further, possibly 
with even higher resolution spectrometers. 
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Fig. A.l. Observed Br y profiles. The observed fluxes have been normalized to the continuum and shifted for easier display. The 
stars are roughly order according to the observed W(Br r ) ofe , from strong emission (left-top) to strong absorption (right-bottom). 
The strongest peak emission (HD 97048) is about 44% of the continuum; the deepest absorption (HR 4796) is 15% of the 
continuum. The four early B stars are displayed in the fifth panel to the right, with a different vertical scale. 



